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Abstract 

This work presents a continuation of our last paper, concerning the theory of the response of an 
antiparallel domain structure in a plate- like electroded sample to external electric field. The theory 
is based on the exact formula for free energy of the system, formed of a central ferroelectric part, 
isolated from electrodes (with a defined potential difference) by a surface layers. Our calculations 
are applicable also to thin films. It is usual to use the term 'extrinsic' for the contribution of 
domain walls displacement to macroscopic properties of a sample. In our last paper we discussed 
the extrinsic contribution to permittivity. In this work we concentrate on extrinsic contribution to 
piezoelectric coefficients in ferroelectrics which are simultaneously ferroelastics. As an example, we 
calculate the extrinsic contribution to piezoelectric coefficient of RbH^PC^, that was recently 
measured in a wide range of temperature below Curie point. 
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I. INTRODUCTION 



Samples of ferroelectric single crystals often posses a surface layers. Its existence greatly 
influences properties of bulk samples 1, 2, 3, 4] as well as of thin films. [5] Equilibrium domain 
structure in the system, mentioned above in the abstract, and the role of the surface layers 
was first discussed by Bjorkstam and Oettel[6j] in a special case of shorted electrodes. In our 
recent paper [7] we reconsidered this problem in a general case of nonzero voltage between 
electrodes, discussing the response of the domain structure to external electric field. Our 
calculations are valid also for thin films and present, in fact, continuation of our discussion 
of domain structures of thin films 1 81. 

n 

In [7] we used our theoretical results for prediction of extrinsic contribution to permittivity 
of the sample. In the next two sections we give a short review of notation, description of the 
model and basic results from[7|. In last two sections we discuss as an example the extrinsic 
contribution to piezoelectric coefficient of RbH^PC^. We compare our predictions with 
the recent measurements of record values ~ 4000 pC/N in temperature range 35 K under 
critical temperature 146 K 9] (see alsojlo|). 



II. DESCRIPTION OF THE MODEL 



We consider a plate-like electroded sample of infinite area with major surfaces perpendic- 
ular to the ferroelectric axis z. Central ferroelectric part with antiparallel domains (2.) is 
separated from the electrodes (0.), (4.) by nonferroelectric layers (l.),(3.) (see Fig. [I]). The 
spatial distribution of the electric field E is determined by the applied potential difference 
and by the bound charge divPo on the boundary of ferroelectric mate- 
rial, where Po stands for spontaneous polarization. Geometrical, electrical and material 
parameters of the system are shown in Fig.fTJ 

We further introduce the symbols 

c = \fe a /e c , g = ^Je a e c , 
and several geometrical parameters: the slab factor 
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2a = a + + a_, 



the domain pattern factor 
2a 

and the asymmetry factor 

a_i_ — a_ 
A = — . 



electrode (z = f/2) 



a+ + a_ 

The ferroelectric material itself is ap- 
proximated by the equation of state 
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where Pq is the spontaneous polarization 
along the ferroelectric axis. This linear 

FIG. 1: Geometry of the model 
approximation limits the validity of our calculations to the temperature region not very 

close below the transition temperature T c . Domain walls are assumed to have surface energy 

density a w and zero thickness. 



III. GIBBS ELECTRIC ENERGY OF THE SYSTEM, EQUILIBRIUM DOMAIN 
STRUCTURE 

Rather cumbersome calculations lead to the following formula for Gibbs electric en- 
ergy per unit area of the system (in Jm" 2 ), which includes the domain wall energy, the 
electrostatic energy whose density is (l/2)i? • (D — Pq) and the work performed by external 
electric sources — VQ, where Q is the charge on positive electrode. 

G = -a R\PA PoA Bt - 2£o£dV (1 + B) 4P ° 2t V Sin2 (wr(1 + A),2) 

7r°™ 2e (1 + B) {Ed + B e c ) e it 2 R{1 + B) ^ n 3 (cothnBR + c/coth nRc) ' 

(1) 

The first term represents domain wall contribution while the last one is the depolarization 
energy. In the second term we recognize the effect of layers (1.) and (3.) and of the applied 
voltage. 

In this model we neglect the mechanical interactions between components of the system. 
For given slab factor B and voltage V, the equilibrium domain structure, characterized by 
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R eq (V) and A eq (V), corresponds to local minimum of G. In general a minimum can be 
found by numerical methods, but for BR ^> 1 and Rc ^> 1, the R eq (V) and A eq (V) can be 
approximated by explicit formulae. For purpose of this paper we use the following formula 
for A cq (V) 

(2) 
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where 



X 



2 \n2{l + Be c /e d ) 



is considered as a small correction and R® q is equilibrium value of R for zero voltage V. For 



the extrinsic contribution to permittivity of the sample we get from (j2J) (see also[7|) 
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IV. EXTRINSIC PIEZOELECTRICITY 



As an example, in this section we work out the approximate prediction of extrinsic con- 
tribution to g?36 of RbH2P04 (RDP), based on our simple model. RDP is a ferroelastic with 
spontaneous strain eo,6, opposite in opposite polarized domains. In the FigJ2h there is x — y 
cut through symmetric domain structure (A = 0, E ext = V/t = 0). The situation after 
application of E ext is shown in Figf2b (i / 0). A simple geometric consideration leads to 
the formula for average strain e e ff,6 of the sample (we neglect the mechanical coupling of the 
central part with the rest of the sample) 



e e ff,6 — £o,gA 



(4) 



For the extrinsic coefficient dze we get from (j2J) and 

-l 
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neglecting the small X and for B < 1 we get more simple formula [l3| 
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V. DISCUSSION 

n n 

An extremely high d 36 under T c for RDP was first reported in[ll|]. In a recent paper [121]. 
Sidorkin deduced the dispersion law of domain wall vibrations, however, in his treatment the 
existence of a surface layers is not explicitly considered. We can fit our theoretical results 
to measured onesfsl - G?36,ext = 4 • 10~ 9 C/N, = 2000 in the 35 K range plateau under 



T c . Using values that roughly apply to RDP: [14J] in (J6J) resp. fl3]) we come to an agreement 
for reasonable value of B = 0.04. Naturally for lower temperatures, the motion of the walls 
is limited by "freezing" and both G? 36iex t and e^ ey ± decrease to zero. It is also interesting, 
for measurements in[9(] with alternating E cxt = 20V/m, that corresponding amplitude of 
alternating A from (j2J) is only 10~ 5 and displacement of the walls with a = 10 /zm is of the 
order 10 _10 m. 



5 



Acknowledgments 



This work has been supported by the Ministry of Education of the Czech Republic grants 
CEZ: Jll/98:242200002 and VS 96006. 



[1] R. C. Miller and A. Savage, J. Appl. Phys. 32 (1961). 

[2] M. E. Drougard and R. Landauer, J. Appl. Phys 30 (1959). 

[3] H. E. Miiser, W. Kuhn and J. Albers, Phys. Stat. Sol. (a) 49 (1978). 

[4] D. R. Callaby, J. Appl. Phys. 36 (1965). 

[5] A. K. Tagantsev, C. Pawlaczyk, K. Brooks and N. Setter, Integrated ferroelectrics 4 (1994). 

[6] J. L. Bjorkstam and R. E. Oettel, Phys. Rev. 159 (1967). 

[7] A. Kopal, P. Mokry, J. Fousek, T. Bahnik, to appear in Ferroelectrics. 

[8] A. Kopal, T. Bahnik and J. Fousek, Ferroelectrics 202 (1997). 

[9] M. Stula, J. Fousek, H. Kabelka, M. Fally and H.Warhanek, J. Kor. Phys. Soc. (Proc. Suppl.) 
32 (1998). 

[10] E. Nakamura, Ferroelectrics 135 (1992). 

[11] L. A. Shuvalov, I. S. Zheludev, A. V. Mnatskanyan and Ts.-Zh. Ludupov, I. Fiala, Bull. Acad. 

Sci. USSR, Phys. Ser. 31 (1967). 
[12] A. S. Sidorkin, J. Appl. Phys. 83 (1998). 

[13] The limit B —* in ([6]) is not correct, because the assumptions needed for validity of ([2]) are 

violated if B is very small. 
[14] P = 4 • 10- 2 C/m 2 , e d = e c = 100, e 0>6 = 0.015, seefl. 



6 



